Localizing the Angular Momentum of Linear Gravity 
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In a previous article [1] we derived an energy-momentum tensor for linear gravity that exhibited 
positive energy-density and causal energy-flux. Here we extend this framework by localizing the 
angular momentum of the linearized gravitational fleld, deriving a gravitational spin tensor which 
possesses similarly desirable properties. By examining the local exchange of angular momentum 
(between matter and gravity) we find that gravitational intrinsic spin is localized, separately from 
"orbital" angular momentum, in terms of a gravitational spin tensor. This spin tensor is then 
uniquely determined by requiring that it obey two simple physically-motivated algebraic conditions. 
Firstly, the spin of an arbitrary (harmonic-gauge) gravitational plane-wave is required to flow in the 
direction of propagation of the wave. Secondly, the spin tensor of any transverse-traceless gravita- 
tional fleld is required to be traceless. (The second condition ensures that local field redefinitions 
suffice to cast our gravitational energy-momentum tensor and spin tensor as sources of gravity in 
a quadratic approximation to general relativity, and simultaneously rids the gravitational field of 
infinite pressure gradients). Additionally, the following properties arise in the spin tensor sponta- 
neously: all transverse-traceless fields have purely spatial spin, and any field generated by a static 
distribution of matter will carry no spin at all. Following the structure of our previous paper, 
we then examine the (spatial) angular momentum exchanged between the gravitational field and 
an infinitesimal detector, and develop a microaveraging procedure that renders the process gauge- 
invariant. The exchange of non-spatial angular momentum (i.e. moment-of-energy) is also analysed, 
leading us to conclude that a gravitational wave can displace the centre-of-mass of the detector; this 
conclusion is also confirmed by a "first principles" treatment of the system. Finally, we discuss the 
spin carried by a gravitational plane-wave. 

PACS numbers: 04.20.Cv, 04.30.-w 
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I. INTRODUCTION 

We recently developed a local description of energy 
and momentum in linear gravity, deriving a gravita- 
tional energy-momentum tensor Tab that describes posi- 
tive energy- density and causal energy-flux [1].^ The pur- 
pose of this present article is to complete our picture of 
local linear gravitational energetics, extending our frame- 
work to quantify the angular momentum carried by the 
field. This approach will localize both the "orbital" an- 
gular momentum and the intrinsic spin of linear gravity, 
the former in terms of Tab, and the latter in terms of a 
gravitational spin tensor s^j,^. Not only is this spin ten- 
sor vital if one is to account for the angular momentum 
possessed by gravity and exchanged locally with matter, 
the formula we derive for it will display a number physi- 
cally desirable algebraic properties, closely analogous to 
those of Tab- 



Armed with a local description of the energy, momen- 
tum, and angular momentum of linear gravity, we will 
be ready to tackle the task of our next paper [2]: to 
understand Tab and 5°^^^ in terms of the familiar theo- 
retical apparatus that has been used to define gravita- 
tional energy-momentum in the past [3-5], and energy- 
momentum in general [6-8]. These developments will 
crystallize the tensors' physical interpretation, deepen 
our understanding of their theoretical underpinnings, and 
suggest a route by which our work might be generalized 
beyond the linear approximation. 

Let us begin by summarizing the key points of the 
program developed in [1].'^ We define the gravitational 
field hab on a flat background spacetimc (A4,gab) by a 
diffeomorphism (p : M ^ Ai that maps the physical 
spacetime {M,gab) onto the background:'^ 



gab — gab + hab- 



(1) 



* l.butcher@mrao.cam.ac.uk 

^ To clarify: we have not performed the demonstrably impossible 
feat of finding a tensor t^^i,, quadratic in the gravitational field 
hab E'nd second-order in derivatives, that is invariant under the 
linearized gauge transformation Shab = ^(a^b)- Rather, we rely 
on a gauge-fixing programme (motivated by key properties of r^t 
and the energetics of an infinitesimal gravitational detector) to 
remove the freedom to perform such transformations, and hence 
arrive at a physically unambiguous description. 



^ As before, we work in units where c = 1, write re = SttG, and 
use the sign conventions of Wald [9]: r;^,^ = diag(— 1, 1, 1, 1), 
[VcVd]';'^ = 2V[,Vd,»;'^ = R^t^y, and R^b = R\^^. We use 
Roman letters (except i,j, k, I) as abstract tensor indices [9, chap. 
2.4] and Greek letters as numerical indices running from to 3. 
The indices i,j, k, I are reserved for spatial components, and run 
from 1 to 3. 

^ As usual, fields defined on M have their indices raised and low- 
ered with Qab, and those on At with g^b- Lorentzian coordinates 
{x'^} are commonly deployed in M, for which = r;^^. 
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The physical spacetimc is assumed to be "nearly flat", 
and (/) chosen such that hab is small everywhere, so that 
the linear approximation to the Einstein field equations 
is valid: 

G,t,''''hcd^Kfab + 0{h^), (2) 

where Tab = 4'*Tab ^ 0{h) is the matter energy- 
momentum tensor Tab mapped onto the background, and 

Gab^'^hcd = ^c^{aK{ - iV^/lab - iVaVf,/! 

+ \9ab {y^h - VcVdh'^) (3) 

is the linearized Einstein tensor G^^fj . 

The gravitational energy-momentum tensor Tab is de- 
fined by seeking a symmetric tensor, quadratic in Vchab, 
which solves 

Vaj/+0*(V,J/)=O, (4) 

neglecting terms 0{h^). In the above equation, J^^"" = 
T\e^ are the (1 energy, 3 momentum) current-densities 
of matter, associated with the (1 timelike, 3 space- 
like) vector fields e^" = {(t)^^)*e^°' , the images of the 
Lorentzian coordinate basis e^" = {d/dx^)'^ that gen- 
erate the translational symmetries of the background; 
the j^"" = T\e^ = t"-^ constitute the energy-momentum 
current-densities of the gravitational field. Consequently 
(4) indicates that the extent to which material energy- 
momentum fails to be conserved at a point in the phys- 
ical spacetime is exactly equal and opposite to the ex- 
tent to which gravitational energy-momentum fails to be 
conserved at the corresponding point in the background. 
Interactions between matter and gravity can then be un- 
derstood in terms of a local exchange of energy and mo- 
mentum between the two. 

It is not possible to construct a Tab to solve (4) for all 
gravitational fields, so a condition must be placed on hab 
in order to proceed. Of all possible symmetric tensors 
Tab, quadratic in S/chab, a-nd all (non-trivial, linear and 
Lorentz invariant) field conditions, only one combination 
solves (4): 

Ufab = jVahcd^bh"'^, (5) 

V^'hab = 0, (6) 

where the overbars signify trace-reversal. Because (6) is 
simply the equation of harmonic gauge, which can always 
be satisfied through a choice of 0, the field condition docs 
not restrict the physical applicability of our approach in 
any respect. In fact, the only effect of the field condition 
is to vastly reduce the gauge freedom in our description 
of gravitational energy-momentum (5). What at first ap- 
peared as a weakness, is in fact a great strength of our 
approach. Essentially, (6) indicates that (j) is to be cho- 
sen such that it maps Lorentzian coordinates {x^} of the 



background onto harmonic coordinates y^{p) = x^{(l){p)) 
of the physical spacetime. This ensures that the energy- 
momentum currents are defined by the generators 
of a harmonic coordinate system; these represent the ap- 
proximate translational symmetries of the physical space- 
time (present due to its small curvature) and give a sen- 
sible replacement for killing vectors in the absence of an 
exact symmetry. 

The gravitational energy-momentum tensor Tab bas 
two notable mathematical properties, in addition to solv- 
ing (4). Firstly, the energy- momentum tensor for any 
(harmonic gauge) gravitational plane-wave 

hab^hab{x''K), k^hab^O, fc^fc,, = 0, (7) 

is completely invariant under the remaining gauge free- 
dom consistent with (6) and (7). Secondly, and most 
remarkably of all. Tab displays the following positivity 
property: all transverse-traceless (tt) gravitational fields 
have positive energy- density and causal energy-flux, for 
all observers. To state this rigorously: if, at some point 
p £ M, the gravitational field hab obeys the transverse- 
traceless conditions 

V'^/la6 = 0, h = 0, U''hab = 0, (8) 

for some timelike vector then Tab satisfies the follow- 
ing inequalities 

v'^Tabv'' > 0, (9) 
V^TacT^v' < 0, (10) 

at p, for any timelike vector u° . 

In order to deal with the last trace of gauge freedom 
that remains after enforcing (6), we examined the energy- 
momentum transferred between the gravitational field 
and an infinitesimal probe, i.e. a matter "point-source" 
with energy-momentum tensor 

Too = MS{x) + \hjd,dj5{x), 

%, = \{hj-L,,)d,5{x), (11) 

Tij ^lijSi^x^ ^ 

derived by shrinking a compact source down to a point.* 
The exchange is rendered gauge-invariant (under the free- 
dom that remains after (6) has been enforced) by the 
monopole-free microaverage: the incoming wave is split 



* M, lij and Lij are the mass, moment of inertia, and angular 
momentum of the source, respectively. Overdots indicate differ- 
entiation with respect tot = x^, and the three spatial coordinates 
are abbreviated x = (x^ , x'^ , x^). In [1], the angular momentum 
of matter was written — Jij ; our change in notation corrects for 
the unusual sign convention chosen in equation (AlO) of [1] , and 
avoids confusion with the energy-momentum current-densities 
7 " 
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into a sum of Heaviside step-functions, and the energy- 
momentum delivered by each is integrated over a vanish- 
ingly smah 4- volume centered on the probe. ^ The result 

{d^T^.)f = -l6ix)h,d,hJJ (12) 

is equal to the bare (i.e. not microavcraged) energy- 
momentum delivered by the incident field in TT-gauge. 
This motivated the program of fixing the final piece 
of gauge freedom by insisting that the incident hab 
be transverse-traceless; consequently, Tab represents the 
gauge-invariant gravitational energy-momentum that is 
accessible to an infinitesimal probe at rest in the tt 
frame. Furthermore, due to the positivity property of 
Tab, this program ensures that the gravitational field is 
always described with positive energy-density and causal 
energy flux. 

The approach we will take for localizing gravitational 
angular momentum will be very similar to the one we 
have just described. Section II of this paper begins with 
the counterpart of (4) for the local exchange of angular 
momentum between matter and gravity. We will show 
that the local change in the angular momentum of mat- 
ter is not entirely accounted for by the change in orbital 
angular momentum 2x[^ry]° carried by the gravitational 
field: gravity's intrinsic spin s°^^ must be included to 
balance the exchange. This argument defines s"^^ up to 
the addition of total divergences, so further requirements 
must be placed on the tensor before we have a unique for- 
mula localizing gravitational intrinsic spin. We achieve 
this in section III by demanding that s^^^, satisfy two 
simple, physically motivated, algebraic conditions, anal- 
ogous to the algebraic properties of Tab- As a result, a 
formula (32) is derived for the spin tensor of the gravita- 
tional field. The gauge freedom of s°^^ is automatically 
nullified by the tt program motivated in [1]; however, 
it is still enlightening to reprise our analysis of the in- 
finitesimal probe and develop a microaverage procedure 
that renders the transfer of angular momentum gauge- 
invariant (within harmonic gauge) without the need to 
fix the gauge completely. This is covered in section IV. In 
section V we examine the role of the non-spatial compo- 
nents of gravitational angular momentum, demonstrating 
that the exchange of non-spatial spin s°gj can displace 
the centre-of-mass of a gravitational probe. We conclude 
our investigation with a calculation and analysis of the 
intrinsic spin carried by a gravitational plane-wave. 

II. LOCAL ANGULAR MOMENTUM 
EXCHANGE 

The purpose of this section, and the one that follows 
it, is to extend the basic framework of [1] to include 



^ Details are to be found in Sections IV B, and IV C of [1]. 



a local description of gravitational angular momentum. 
Unlike our work on Tat-, eliminating gauge freedom will 
not be a major concern: we already know that harmonic 
gauge (6) is necessary (thus we shall enforce this condi- 
tion throughout the article), and that the last trace of 
freedom must be removed by insisting that the incident 
field be transverse-traceless. We begin by formulating 
the local exchange of angular momentum. 

As noted in Section I, the material energy-momentum 
current-densities J^" are formed by contracting T^^ with 
the vectors e^'' = {(f>~^)*e^^, the pushforwards of which 
under cf) generate the translational symmetries of the 
background. Therefore, to define material angular mo- 
mentum current-densities J^y"', we must contract T\ 
with the vector fields (0~^)*(2a;[^ei,]''), the pushforwards 
of which generate the rotational symmetries of the back- 
ground:^ 

j^/ ^ n(r')*(2x[^e,]'') 

= 2ny[^e,]''. (13) 

As usual, {x^} comprise a Lorentzian coordinate system 
on the background, and y^{p) = x^{4){p)) are the image 
of these coordinates in the physical spacetime. The {j/^} 
are harmonic (that is, V^y'^ = 0) as a result of the gauge 
condition (6). 

We wish to explain the effect of the gravitational field 
on the angular momentum of matter in terms of a local 
exchange of angular momentum between the two. Just 
as (4) captured this idea for energy-momentum, we will 
require 

Va.7^."+0*(VaJ^.")=O (14) 

for angular momentum, where j^j^" is the angular mo- 
mentum current-density of the gravitational field. Ne- 
glecting terms 0{h^), equation (14) is equivalent to 

VaJ^/ = -0*(Va(V)) 

= -0*(nVa(2y[^e,]'')) 

+ 2Va(x[^e,]'')] 
= -n(V,/i,'')x[^e,]- 

-2(n-/.-fe6)Va(x[^e,]''), (15) 

where in the last line we used (j)*{T\) = <!)* [g"-''T^b) = 
f% - h'^'fcb + 0(/i^) and fab = fba- As we now have an 



® We use the term "rotational symmetry" here as a shorthand for 
both rotations and Lorentz boosts. The three independent vec- 
tor fields 2a;[iej]" generate rotations (so that 23;[ie2]'' rotates 
about the x^-axis, for example) and hence define angular mo- 
mentum current-densities. The 2a;[gej]'' generate boosts (in the 
x'-direction) and define moment-of-energy current-densities, the 
interpretation of which we explore in section V. 
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equation relating tensors defined on the background, we 
can express these tensors in terms of their components 
in the Lorentzian coordinate system: 

^-x^^f'-f'd,^K0 + 2hp^;r,f. (16) 

Finally, we recall the field equations (2) in harmonic 
gauge, 

d^Kb = -2kT„^, (17) 
and eliminate Ta^ from (16): 

^d^ [ix^^T^^^ + hp^.d^h/ln] . (18) 

This is rather surprising result, and one that reveals the 
importance of the gravitational field's intrinsic spin. The 
first term in the square brackets clearly represents the or- 
bital angular momentum of the field: it takes the familiar 
form X X p and is the result of the tangential linear mo- 
mentum about the origin. The second term, in contrast, 
does not depend explicitly on x'*; it measures the extent 
to which the field itself is spinning at a particular point, 
and contributes the same gravitational angular momen- 
tum without regard to where this spin is taking place. 
We are forced by (18) to accept that the angular mo- 
mentum of the gravitational field is not simply orbital, 
but also has an intrinsic component: 

where s^^j^ is the gravitational spin tensor (composed of 
intrinsic spin current-densities) without which the local 
exchange of angular momentum would not balance. Of 
course, the division of angular momentum into orbital 
and intrinsic components is not a new idea, and the form 
of equation (19) originates from standard flat-space field 
theory [10, 11]. In general, the Noether current of a ro- 
tational symmetry cannot be constructed entirely from 
Noether currents of translational symmetries: the mis- 
match, born of the field's tensorial (or spinorial) struc- 
ture, is called intrinsic spin.'' More neatly, and of greater 
relevance to our later analysis, the energy-momentum 
tensor and the spin tensor can be derived separately from 
a Lagrangian by 'gauging' the translational and rota- 
tional symmetries of spacetime and taking the functional 



Essentially, this is because a tensor field undergoes two types 
of transformation when it is rotated. A vector field A^^{x), for 
example, becomes A^^A''{A~^{x)); in the parlance of quantum 
field theory, this can be understood as a displacement x A{x) 
generated by the orbital angular momentum operator x X p, and 
a pointwise Lorentz transformation A'^ — > A^^A'^ generated by 
the spin operator. 



derivatives with respect to the two gauge fields. In the 
paper that follows [2] we construct r^,y and s"^^ accord- 
ing to this method, confirming that our formulae for r^^ 
and s"^^ (soon to be derived) are in keeping with the 
established concepts of energy- momentum and spin.® 

III. GRAVITATIONAL INTRINSIC SPIN 
TENSOR 

Our immediate goal, of course, is to arrive at a formula 
for s"^^^ in terms of hap. With this in mind, it is tempting 
to solve (18) simply by setting 

Ks'^^, I hp^^d'^h/, (20) 

and declare that we have found our local description of 
gravitational spin. However, this is not the only solution: 
the exchange equation (18) only defines s"^^ up to terms 
with identically vanishing divergence, so further demands 
must be made of the spin tensor before it can be deter- 
mined uniquely. Obviously, s"^j, should have the same 
basic properties as r^i/: it should be a local, quadratic, 
Lorentz-covariant function of hap , and contain no dimen- 
sionful constants other than k.^ The general solution to 
(18) is then 

Ks"^, = h0[,d''h/ + dpj:''P^,, (21) 

where T,"^ ^i, is any local, quadratic, Lorentz-covariant 
function of hap (but not its derivatives) that obeys 

^ — ^ /if — ^ f/j- \^^) 



' Because spin tensors are usually associated with asymmetric 
energy-momentum tensors, it is worth mentioning that the sym- 
metry of Tfj^y does not contradict the existence of s'^^^j. Typ- 
ically, one argues that Tf^^^j ^ describes finite torques act- 
ing on infinitessimal regions [12], and then states that this is 
only acceptable if one can interpret these torques as generating 
intrinsic spin: das""^^ = 2t[^^] . Clearly, this argument does 
not run in reverse: the presence of a spin tensor does not re- 
quire that the energy-momentum tensor be asymmetric. A sym- 
metric gravitational energy-momentum tensor simply indicates 
that there are no torques on infinitesimal regions due to gravity, 
and so (in the absence of matter) the spin-tensor is conserved: 
daS"^^ = 2t[^^] = 0. 

^ This last stipulation (which forces the terms in s",,^ to contain 
exactly one derivative, in order that they have the correct units) 
is essentially unavoidable within the context of classical general 
relativity: k is the only dimensionful constant available. If we 
allow ourselves to use Planck's constant h (as we would for a 
quantum theory) or introduce a new dimensionful gravitational 
constant (as would arise in a higher-derivative theory of gravity) 
then higher derivative terms would be dimensionally permissible 
within the spin tensor; nonetheless, these higher-derivative terms 
would each be multiplied by small factors (such as the Planck 
length) that would ensure the terms were negligible within the 
low-curvature regime of the theory that corresponds to classical 
general relativity. 
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The most general tensor that can be formed from hap 
this way is 

+ <5p^<5f] {A^h^ + A^h^sh^') , (23) 

where the are arbitrary dimensionless constants. 

Equations (21) and (23) describe the range of possible 
gravitational spin tensors that account for the angular 
momentum exchanged with matter; the aim of this cur- 
rent section is to find a distinguished member of this set, 
deserving of its physical interpretation. 

We encountered a similar "superpotential" ^'^ freedom 
when deriving r^^ in [1] , and extinguished it immediately 
by insisting that the energy-momentum tensor should be 
free of second derivatives. Unfortunately, this tactic is of 
no use here: all the terms in s"'^^ have the same form 
hdh, and so cannot be distinguished from one another 
by their differential structure. Instead, we must place al- 
gebraic requirements on the spin tensor, and we shall do 
so by choosing two conditions that are physically well- 
motivated, and closely analogous to the algebraic prop- 
erties of T^i,. 

A. The Plane-wave Condition 

Condition 1: The spin tensor of any (harmonic gauge) 
gravitational plane- wave (7), with wave- vector A;'^, must 
obey 

(X fc". (24) 

Clearly, this ensures that spin flows in the direction of 
propagation of the wave, a physically reasonable request 
that reciprocates the property t^i, oc k^k^, of plane- wave 
energy-momentum. Substituting (7) into equations (21) 
and (23), we find that the condition (24) holds for all 
harmonic gauge plane-waves if and only if 

A2 = -Ai, Ai = A^/A, A3 = A = 0. (25) 

This leaves us with a much smaller range of spin tensors 

^s-^^u = /^m.a"/!^]^ -t- A^dp {h"[^h^.]) , (26) 

parametrized by ^i. 

Of course, the influence of this first condition is not 
limited to gravitational plane- waves. In fact, the restric- 
tion (26) automatically endows the spin tensor with two 
highly desirable properties that apply to much more gen- 
eral gravitational fields. Furthermore, one can check that 
these two properties occur only if the spin tensor takes 



These superpotentials arc so ealled because they are total deriva- 
tives. They bear no relation to the homonymous concept from 
supersymmetric field theory. 



the form (26); hence the logic can be reversed, with both 
properties taken together as conditions on s"^^, and (24) 
derived as a consequence. 

Property la: The spin carried by a transverse-traceless 
gravitational field (8) is purely spatial: 

/lOc =0, h^O, d^h.j =0 s"o, = 0. (27) 

Not only are transverse-traceless fields blessed with pos- 
itive energy-density and causal energy-flux, now we see 
they carry only standard spatial spin! This result is akin 
to the Frenkel condition [13] that constrains the spin- 
tensor of a Weyssenhofi^ fluid [14, 15]: S^q^ = 0, in the 
rest-frame of the fluid. The only difl[erence here is that 
the gravitational fleld, being masslcss, has no rest-frame; 
in its place, the TT-frame defines the space/time split. 

The reader should not be under the impression that the 
non-spatial spins s^q, are completely unphysical, how- 
ever; as a matter of fact, they have a simple physical 
interpretation. In section V, we explain that the non- 
spatial angular momentum current-densities jg," local- 
ize gravity's Moment-of-Energy, the conserved quantity 
associated with the symmetry of the background un- 
der Lorentz boosts. Accordingly, the intrinsic current- 
densities s^Qj signify an "internal displacement of en- 
ergy" of the field. This alters the gravitational Moment- 
of-Energy just as the "internal spinning motion" , signi- 
fied by s"jj, contributes to the total gravitational an- 
gular momentum. Due to (27) it is now clear that the 
transverse-traceless field does not carry these internal dis- 
placements, and hence, that the location of gravitational 
energy is determined by r^i, alone. The s"gj still play 
an important role in the local exchange of Momcnt-of- 
Energy with matter (see §VB) because TT-gauge cannot 
be adopted where T)ii/ ^ 0. 

Property lb: All static distributions of matter give rise 
to spinless gravitational fields: 

f^, = p{x)SX => sV=0- (28) 

The meaning of this statement is intuitively obvious: 
matter must be in motion if it is to generate gravita- 
tional intrinsic spin. It is worth remembering that the 
linearized gravitational field due to static matter is just 
the Newtonian potential $, so (28) is equivalent to the 
statement that $ has no spin. This is in keeping with 
our observation in [1] that the gravitational field corre- 
sponding to a static Newtonian potential has the energy- 
momentum tensor of a masslcss scalar field. 



The WeyssenhofT fluid is simply a perfect fluid writh intrinsic spin. 
Note that the massive spin-1/2 field (described by the classical 
Dirac Lagrangian) also obeys the Frenkel condition, if one takes 
the charge current-density to define the field's 4- velocity [11]. 
In order that the distribution does not collapse under its own 
gravity, the matter will also have stresses Tij ~ 0{pK) ~ 0{h?), 
but these can be neglected in the linear approximation. 
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B. The Traceless Condition 

So far we have placed one algebraic condition on the 
gravitational spin tensor and removed all but one of the 
superpotential degrees of freedom. Our second condition 
will fix Ai and determine s"^^^ uniquely. 

We begin by noting that (as shown in our subse- 
quent paper [2]) the superpotential 9,3(/i"[^/i'',y]) asso- 
ciated with Ai plays a distinguished role in the quadratic 
approximation to the Einstein field equations. To avoid 
a major diversion, let us make the following statement 
here and postpone its proof until [2]: if the difference 
between the physical metric and the background metric 
4'*9ab ^ gab is a local function of the gravitational field 
/lab, then the quadratic approximation to the vacuum 
Einstein field equations can be written as 

G^^'-^Kp = [v. + + s.^" - s"^J/2] , (29) 

if and only if 

Ai = -1. (30) 

The tensor in square brackets on the right-hand side of 
(29) is the Belinfante energy-momentum tensor of the 
gravitational field, combining t^i, and s"^^, into a single 
object. According to (29) this tensor acts as a source 
for the gravitational field, fulfilling the role played by 
the (Belinfante) energy-momentum tensor of matter T^i, 
at linear order in the no n- vacuum equations (2). If we 
wish to be able to interpret r^^jy and s"^^ as a genuine 
localisation of gravitational energy-momentum and spin, 
then not only must they (i) account for the exchange of 
energy-momentum and angular momentum with matter, 
but they must also (ii) generate gravity alongside matter 
in the Einstein field equations; for this reason we must 
insist that the remaining superpotential freedom be ex- 
tinguished by setting Ai = —1. Were any other value 
of Ai to be chosen, then a non-loeal field redefinition 
hi^u — > h^i, -\- 0{h?) would be needed to bring the field 
equations into the form (29), and the physical metric 
would no longer be a local function of h^i,. 

Setting Ai = —1 also has an important effect on the 
spin tensor (26) independent of its role in the quadratic 
field equations. Once the gauge has been fixed, and is 
transverse-traceless, Ai — ~1 guarantees that the trace 
of the spin tensor will vanish: s"^,^, = 0. In fact, because 
the spin tensor (26) has this property only if Ai = —1, 
it is possible to fix the final piece of superpotential free- 
dom by placing a second algebraic condition on the spin 
tensor, as follows. 

Condition 2: The spin tensor of a transverse-traceless 
gravitational field (8) must be traceless: 

hoa. =0, /i = 0, ^^h,J =0 0. (31) 

In appendix B, we further explore the physical interpre- 
tation of this condition. As we argue therein, the trace 
of the spin tensor is proportional to an infinite pressure 



gradient carried by the gravitational field. Thus, on the 
understanding that such pressure gradients are unphys- 
ical and must be avoided, this discussion offers an in- 
teresting alternative justification for the traceless condi- 
tion. It is also worth remarking that the above condition 
strengthens the similarity between gravitation spin and 
standard examples of material spin: the spin tensors of 
the Weysscnhoff fluid [14, 15], and the spin-1/2 field [11], 
are also traceless. 

By design, the spin tensor (26) is consistent with Con- 
dition 2 if and only if = — 1; as a result, we arrive at 
our final formula for the gravitational spin tensor: 

Ks^^,^2hp[,d^"h^]Pl (32) 

This is the unique local, quadratic, Lorentz-covariant 
function of /i^j^ that accounts for the local exchange of an- 
gular momentum with matter (14) in harmonic gauge (6), 
satisfies the two physically well-motivated algebraic con- 
ditions (24) and (31), and contains no dimensionful con- 
stants other than k}^ This is an exceptionally compact 
formula, and one that embodies a remarkably parsimo- 
nious description of gravitational spin: for a transverse- 
traceless field, s"^^ is specified by no more than 9 inde- 
pendent components (due to (27) and (31)) as opposed 
to the 24 that would be needed in the generic case. 

This completes the foundational portion of the article. 
Following the structure of [1], our next task is to apply 
our newly assembled framework to an investigation of the 
angular momentum absorbed by an infinitesimal gravita- 
tional detector. Section IV will focus on the exchange of 
standard (i.e. spatial) angular momentum jij" , and the 
microaverage that renders this process gauge- invariant; 
section V concerns the interpretation of non-spatial an- 
gular momentum jjg", and the physical consequences of 
its exchange. A reader whose primary interests are the 
theoretical underpinnings of r^^ and s"'^^, may wish to 
skip to [2] at this point: knowledge of sections IV, V, and 
VI will not be necessary for the discussion therein. 

IV. ANGULAR MOMENTUM 
MICROAVERAGE 

Having derived the formula (32) for gravitational spin, 
we now possess a complete description of the local en- 
ergy, momentum, and angular momentum carried by the 
linear gravitational field. Our first application of this 
framework will be an analysis of the angular momentum 



This derivation has taken place entirely within harmonic gauge, 
so it goes without saying that s"ji^ has only been uniquely de- 
fined up to the addition of terms proportional to d^^h)^^. As 
long as we remain in harmonic gauge (which we must if we are 
to interpret r^^ and s"^^ physically) then such terms can clearly 
be ignored. For a discussion of the unique extension of t^^ and 
s"^^ beyond harmonic gauge, see section II of our subsequent 
paper [2], 
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exchanged with an infinitesimal probe. This will allow us 
to revisit the microaverage, the procedure which defined 
the gauge- invariant energy-momentum transferred onto 
the probe/* and motivated the (equivalent) program of 
preparing the incident field in transverse-traceless gauge. 
Clearly, this gauge-fixing program also provides us with 
an unambiguous definition of the angular momentum ex- 
changed with the probe. What is not obvious, though, is 
whether a microaveraging procedure can also achieve this 
effect, allowing us to define a gauge- invariant exchange 
of angular momentum that does not rely on gauge-fixing. 
The aim of this section is to confirm the truth of this idea. 

We shall consider a system that is almost identical to 
the one described in section IV A of [1]: a point- like 
detector in the path of a gravitational "pulse" plane- 
wave. The gravitational detector will once again con- 
sist of an infinitesimal point-source at x ~ 0,^^ the 
energy-momentum tensor of which is given by (11) as 



0. The gravitational field 



T 7^ wave I source 



(33) 



is the sum of the incoming gravitational wave. 



7^ wave 



(1,-1,0,0), k^'A 



const., 

0, (34) 



and the field 



generated by the detector. 



g^^souree 



(35) 



It is important to recognise that the plane-wave (34) is 
not quite the same as the one we used when defining 
the energy-momentum microaverage. There, the gravi- 
tational wave had the profile of a Heaviside step func- 
tion H, and this brought about an exchange of energy- 
momentum d^Tfj^u ~ dhd^h ~ S{t — to)S{x) that was 



^'^ As our description of gravitational energetics only exists in har- 
monic gauge, we need only consider gauge transformations which 
do not break the harmonic condition (6). Hence, we use the term 
gauge-invariant to mean invariant under the gauge freedom that 
remains after enforcing the harmonic condition. 
It might appear that we risk a loss of generality in placing 
the probe at the origin, but this is not the case. To ex- 
plain, let us consider a uniform translation of the coordinates 
x'' — > x'' + a^; the probe then lies at a; = a, and according 
to (19) the only effect on the gravitational angular-momentum 
current-density is ^j^j" = ^^[iTj]". Because a* is constant, the 
exchange of angular-momentum associated with this term is sim- 
ply da{i^j^j°') = 2a^fd"Tj-^a, and we already know from (12) that 
d^Tia (which quantifies the local exchange of linear momentum) 
is rendered gauge-invariant by the monopole-free microaverage: 
(9a(Aj^^"))j^^= 2a[j (cJ^Tj] a) j^^ Clearly, this term accounts for 
the angular momentum that results from the transfer of linear 
momentum onto the detector; by assuming that the probe is 
at X = in what follows, we are simply ignoring the trivial ex- 
change of angular momentum associated with the detector's bulk 
motion. 

^® We take this limit as the size of the source shrinks to zero. The 
detector is then a form of generalised "test-particle" with negli- 
gible self-interaction in comparison to the effect of the external 
field. 



confined to an infinitesimal spacetime region over which 
we could average. The same is not true of angular mo- 
mentum, however: a step function wave will give rise 
to a local exchange daj^^" including a term Sas"^^ ^ 
hd^h ~ H{t — to)S{x) that is not localized a.t t ~ to 
and is therefore unsuitable for microaveraging. We have 
no choice but to use a delta-function wave to generate a 
point-like angular momentum exchange. This will be 
the only modification needed to adapt the microaverage 
for angular momentum. 

Following the same reasoning that took us to equation 
(42) of [1], we find that the exchange of spatial angular 
momentum for this system is given by 

= -^k[jXi]S{kax" - to) IkiAkiS{x) 
-2{iki~Lki)AkodiS{x) 
+ {2MS{x)+hidkdiSix))Aoo 

- S{kaX°' - to) Ao[^{iJ]k - Lj]k)dkS{x) 

- Ak[tlj]k5{x) . (36) 

As was the case with energy-momentum, the local ex- 
change of angular momentum (36) is clearly not invariant 
under the gauge transformations 



SA 



E(^k^), 



E„ 



const. 



(37) 



which neither break the harmonic condition (6) nor alter 
the form (34) of the wave. This gauge dependence can 
be dealt with in one of two ways. The simplest approach 
is to invoke the familiar tt program, insisting that the 
incident field be transverse-traceless: = A^'J). The 
alternative, which we will now examine, is to integrate 
over the infinitesimal interaction region and render the 
exchange gauge-invariant without gauge-fixing. The two 
methods give identical results, as we shall soon show. 

The microaverage (. . .)ta is defined, just as it was in 
[1], by 



{f)t,=5{x)5{t^to)Yiu, 



fd^x, 



8. (to) 



where Be{to) = {(t, x) : \t - to\ < e, \x\ < e}. (38) 



^'^ One might try to use a pulse based on derivatives of the delta- 
function, but the process of splitting a general wave into such 
pulses is non-local and introduces an arbitrary constant of inte- 
gration. 

The lesson here is that the microaverage is not a process in which 
we split the incident wave into a particular sort of pulse: as we 
have seen, the profile of the pulse depends on what exchange we 
arc microaveraging. Rather, it is a process in which we split the 
wave such that the local exchange (of energy-momentum 9ct r , 
or angular momentum dajij") takes a particular form: a series of 
delta-function pulses (and possibly derivatives of delta-functions) 
each of which can then be averaged over a vanishingiy small 4- 
volume. 
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Applying this definition to (36) and integrating by 
parts/^ we arrive at 

{daj^")tQ 5{x)5{t - to) ky{i,]kAkO + /jJlAoo) 



(39) 



Although it is far from obvious in its current form, this 
equation is in fact invariant under the gauge transforma- 
tions given in equation (37). The easiest way to demon- 
strate this is to examine each component in turn and to 
use k^Afj^^ = in the following form: 



^00 + An + 2Aoi = 0, 
Ao2 + Ai2 ^ 0, 



A22 + A33 = 0, 
Ao3 + Ai3 = 0. 



(40) 



After a great deal of canceling, one finds that 

(9aJ23">to = Six)Sit - to) (^+/23 - A^ (/22 - /33)/2 
(9aJl2">to = -S{£)S{t - to) (A+/i2 + ^x/13) /2, 
(5aji3")to = -Simt - to) (A^h2 - A+/13) /2, (41) 



all of which depend only on the transverse components 
of the wave -Ax 

= ^23, A+ = (^22 - ^33)72 which are 
invariant under (37). Considering that the microaver- 
age was developed purely for the purposes of energy- 
momentum exchange, it is gratifying to discover that 
it renders the exchange of angular momentum gauge- 
invariant as well. 

It is possible to write the above relations (41) in a more 
compact form: 



{daji")to = S{x)S{t - to)AllIj^k, 



(42) 



where is the transverse-traceless part of A^^, the 
only non-zero components of which are A22 = —^33 — 
A+ and A2I = ^32 = ^x- As previously advertised, 
this is exactly the same result as would be obtained from 
applying the tt program to the bare angular momentum 
exchange (36): 



^ + AllSiKx'^ - to)] 

= 5{x)6it-to)AH,i,],. 



(43) 



The only subtlety with this calculation is that one must 
set Xi5{kaX°' — to)S{x) = 0, which is valid as an identity 



between distributions on test functions that are differen- 
tiable with respect to t at {to, 0). 

The angular momentum microaverage need not be re- 
stricted to plane- wave pulses: we can generalise equa- 
tion (42) following the same procedure as the energy- 
momentum case. First we note that an arbitrary 
(harmonic-gauge) plane-wave 



S^,(fc„x"), fc^'B^, =0, (44) 



can be split into a sum of individual pulses 

/oo 
B^,{to)S{k^x^^to)dto, (45) 
-00 

and the angular momentum exchange of each pulse mi- 
croavcraged separately:^^ 



/ ■ CKr 7 source I T wavciX 



+ B^,ito)Sika.x'' ~ to)])todto. 



(46) 



Second we recall that any incident field ft-Jfj, can be ex- 
pressed as a sum of plane- waves, at least locally. Because 
(42) is linear in the incident field, we can split any inci- 
dent field into a sum of plane- waves, each of which can be 
split into a sum of pulses, then perform the microaverage 
on each element and reassemble the result. The general 
formula is therefore 



(5aj,y")/5 = S{x)hllij^k, 



(47) 



where /i^^ is the transverse-traceless part of h'^^^. 

This concludes our analysis of the spatial angular mo- 
mentum transferred onto the probe. The non-spatial cur- 
rents jgj" can also be absorbed by the detector; the ex- 
change equation (14) then ensures that the shift in grav- 
ity's moment-of-energy is accompanied by a displacement 
in the detector's centre-of-mass. This is a rather surpris- 
ing phenomenon, and one that, to our knowledge, has not 
been discussed in the literature. Under resonant con- 
ditions, this effect can cause the detector to "walk" in 
a direction transverse to the gravitational wave.^^ The 
next section is devoted to a detailed examination of this 
phenomenon. 



For each term, integrate by parts to move derivatives from S{x) 
onto the Xi5{kc,x°' —to) or S{kaX°' —to) part of the term, convert 
diS{kaX°' — to) = —5iiS{kax" — to), and integrate by parts once 
again to send the time-derivatives to the M, Jij , lij part of the 
term, recalling that M = Jij = 0. Note that at least one of the 
spatial derivative must act on the Xi in front of the orbital terms: 
those terms where Xi is left untouched will vanish because S{x) 
will set Xi = when the integral is finally evaluated. 



This microaverage carries the subscript fS to remind us that the 
wave has been split into 5-function pulses, rather than Hcaviside 
steps. 

This should not be confused with the motion associated with 
the linear momentum that the probe gains according to (12). 
There, a resonance between the detector and the incident wave 
gives rise to a longitudinal acceleration, and the velocity gained 
in this process remains after the has wave passed. 
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V. MOMENT OF ENERGY 

Through its unification of space and time, and energy 
and momentum, special relativity fused together the once 
disparate notions of angular momentum and centre-of- 
mass. In this section we review this idea in terms of local 
currents, and offer an interpretation for the non-spatial 
intrinsic spin currents s"oi ■ We also examine the local 
exchange of moment-of-energy between the gravitational 
field and an infinitesimal detector. In appendix A we 
confirm that this phenomenon is also predicted by a "first 
principles" description of the system. 



A. Definitions and Interpretation 

It goes without saying that the non-spatial components 
jgj" and Joi" are needed to form the Lorentz-covariant 
currents j^j," and J^^"; thus, at the most basic level, 
these non-spatial components carry the interpretation of 
standard "spatial" angular-momentum as seen by a mov- 
ing observer. Beyond this, the non-spatial components 
carry an additional interpretation that is quite distinct 
from spatial angular momentum. They are the current 
densities of a conserved 3-vector quantity: the Moment- 
of-Energy. 

To explain, let us first define the momcnt-of-cnergy Xi, 
total linear momentum Pi, and total mass/energy M for 
matter: 




P, 



noting that the centre-of-mass x^"'' is simply the moment 
of-cnergy normalised by the total mass/energy: 



(48) 



xf^ = XJM. 



(49) 



The total non-spatial angular momentum of matter is 
then 



— 7 i3 



= -tP,+X,, (50) 
where we have written y'^ — t.^"^ In the absence of the 



Note that we use the same symbol t to represent the value of the 
time coordinate jy" in physical spacetimc and the time coordinate 
of the background. This has the advantage of allowing us 
to drop the distinction between the physical quantities Xi, Pi, 
M, x^^\ and their background representations 4>*{Xi), (p*{Pi), 

<i>*{M), <i>*{xf^): the first set are functions of only, the second 
set of only, and the two sets are numerically equal when x^ = 
yO. 



gravitational field (/i^i/ = 0) the angular momentum cur- 
rents arc conserved, 

da{V^J,,.n^V^^aJ^,/ ^0, (51) 

and as a result, 

dt (X, -tPi)^ 0. (52) 

Furthermore, the conservation of energy-momentum 
(5„(y^T"^) = y^aJ^" = 0) ensures that = 0, 
and leads to the following global conservation law: 



X^-P^ 0. 



(53) 



This equation integrates to Xi ~ tPi + Xi\t=Oi which on 
substitution into (50) gives 



gJo,"d^y = X,\t^o, 



(54) 



which is constant by definition. In other words, the total 
non-spatial angular momentum is equal to the moment- 
of-energy at t = 0, Si conserved quantity which we will 
refer to by the acronym MoE, where the stipulation "at 
t ~ 0" should be taken as given. 

The same analysis can be performed for the gravita- 
tional field in the absence of matter. Working in the 
background, we define 



TQox^d X, 
s\.^d^x. 



r°d\. 



X, = XI + XI 



(55) 



Then the total non-spatial gravitational angular momen- 
tum is given by 



j Jo.Od^a; = X,- tr^ 



(56) 



which, due to dajoi" = and 9a t" = 0, is conserved: 

X, - Vi = 0, (57) 
J jo,"d^x = X,\t=o- (58) 

We conclude from this that the Jq," are the current- 
densities of the conserved quantities Ai|t=o that consti- 
tute the gravitational MoE. 

As (55) makes clear, the non-spatial spin densities s^Qj 
shift the gravitational MoE by Xf, displacing it from 
the value X[ that would have been expected from tqo 
alone. This suggest that the s^q^ represent an "internal 
displacement of energy" at a point (analogous to the no- 
tion of s^jj as "internal spinning motion" at a point) so 
that the field's energy lies locally ofF-centre. The value of 
Too(p) still represents the density of gravitational energy 
at the point p, but an asymmetry in the distribution of 
the energy "within the point", quantified by s'^gj. shifts 
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the MoE by a small amount. Because s°gj = for any 
transverse-traceless gravitational field, these internal dis- 
placements rarely arise when describing the energetics of 
the gravitational field in vacuum. However, as TT-gauge 
cannot be adopted where T^y ^ 0, the s'^Qj inevitably 
play an active role in the exchange of MoE between mat- 
ter and gravity. 



B. Moment of Energy Exchange 



When matter and gravity interact, neither jg^" nor 
Jgj" are independently conserved, and MoE is exchanged 
between them according to (14). Consequently, the con- 
servation laws (53) and (57) are broken. 



X, - = AX, ^ 0, 
X, - = AA-, ^ 0, 



(59) 
(60) 



but the extent to which they are broken is exactly equal 
and opposite 



AA^, + AX, = 0. 



(61) 



To understand this process in general, we turn once 
again to our preferred testing ground; an infinitesimal 
detector in the path of a gravitational plane- wave. Unlike 
our analysis of angular momentum for this system (§IV) 
we will not employ the microaverage here. The reason 
for this is simple: the microaverage does not produce a 
gauge-invariant description of the exchange of MoE. In 
contrast to angular momentum and energy-momentum, 
the gauge-invariant modes of the gravitational field do 
not deliver MoE evenly across the whole detector, they 
are biased by a dipole term proportional to di5{x).^^ The 
microaverage is therefore unable to capture the exchange 
properly, as it can only produce quantities proportional 
to 5{x). This is a notable qualitative difference between 
the exchange of angular momentum and MoE, but in 
reality it poses no practical difficulty: we can still remove 
the gauge dependence by insisting that the incident field 
is transverse-traceless. 

With this in mind, we consider the same system as 
described in section IV with one exception: the incident 
field is an arbitrary transverse-traceless plane- wave. 



BH = B^^=0, (62) 



This pointwise internal structure (spinning motion and displace- 
ments) presumably takes place in the tangent space of the man- 
ifold, where the gravitational field is defined. 
This global exchange equation follows directly from the local 
exchange equations: multiply equation (14) by 0* {\/—g) = 
yj—g + 0{h), discard terms OQv^), and integrate over the spa- 
tial coordinates. This gives AA'i -|- AXi — tij'i + Pi) = 0, and 
7^1 + Pi = follows from the local exchange of linear-momentum 
(4) by exactly the same method. 
This can be seen in equation (64) below. 



rather than a pulse. Taking the same steps that were 
used to derive (42) of [1], and deploying the distributional 
identity Xi6{x) — 0, we find that the local exchange of 
non-spatial angular momentum is 

<9aJio" = ^<9"t„i, (63) 
i = 2, 3 : a„j,o" = B]^ (4 - Lfe,) d,5{x)/2. (64) 

As the longitudinal (63) and transverse (64) equations 
represent two very different phenomena, we shall examine 
them separately. 

Equation (63) is essentially trivial: it accounts for the 
extra MoE that arises from the exchange of linear mo- 
mentum in the direction. To demonstrate this, let us 
take the time derivative of (56): 

X,-V.,-m^ JdoJo.,°d^x = y^^jg^^d^x. (65) 

Unlike the non-interacting case, we now have 

P.- JdoAd^x = y"a„r",d\, (66) 

which is nonzero in general. Consequently, 

AX, = X-r,^ j tda.T", + 9ajo^"d'^. (67) 

Thus, the quantity that describes the local exchange of 
MoE is in fact the sum 

taoT", + a„jo,", (68) 

as it is this combination which contributes the extra in- 
crease in Xi beyond what would be expected from sim- 
ply integrating Vi^t) with respect to time. Because the 
gravitational wave only deposits momentum in the lon- 
gitudinal direction (see equation (12)) this argument has 
no effect on the interpretation of (64); however, equation 
(63) reveals that 

td^T\ + 9„jgi" = td^T\ + {-td^T\) = 0, (69) 

confirming that there is no exchange of MoE in the x^- 
direction, only the exchange of linear momentum. The 
center-of-mass of the detector will accelerate in the x^- 
direction, but this acceleration will be exactly what one 
would expect from the linear momentum transfer dis- 
cussed in [1]. 

In comparison, the exchange of transverse MoE (64) is 
considerably less trivial. The first complication is that 
cJctJio" indicating that the transfer of MoE 

occurs within a dipole-like distribution, taking opposite 
signs at opposite ends of the detector. In general, these 
effects will partially cancel each other, so a more per- 
tinent quantity to calculate (rather than the local ex- 
change) is the total MoE exchange over the whole detec- 
tor: 

Air, = -AA", = - jtdo^T",, + a«ig,"d3x 

= i3-(t)(Ai-L,i)/2, (70) 
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for i = 2,3. This equation describes the transverse drift 
in gravitational MoE, and via (61), the opposite drift in 
the matter MoE. 

In general, the centre-of-mass of the detector (49) will 
move according to 

if ^ = (Air, + P,)/M - X,M/M\ (71) 

under the influence of the gravitational wave. Focus- 
ing our interest on the transverse directions (for which 
Pi = for all time) we note that the last term in (71) is 
the product of two small quantities {Xi and M) and can 
therefore be neglected in comparison to the first term, 
which only contains one small quantity (AX).'^^ Making 
these simplifications, and substituting (70) into (71), we 
finally arrive at a formula for the transverse motion of 
the detector's centre-of-mass: 

z = 2,3: if ^Bj^[h,-Lki)/2M. (72) 

It is important to realise that this motion is not simply 
a "coordinate effect" . If we were to place a free particle 
at rest at the origin, then because the plane-wave is tt, 
this reference point will remain at a; = indefinitely. 
Equation (72) therefore predicts the displacement of the 
centre-of-mass relative to this reference point, and the 
proper distance between the two points will be, to lowest 
order, equal to the Euclidean distance in the background. 

In passing we also note that, when = 0, the accel- 
eration of the centre-of-mass is exactly that of a spinning 
test-particle (of mass M and angular momentum Lij) as 
predicted by the linearized Papapetrou-Dixon equations 
[16, 17] in transverse-traceless gauge: 

Mif^ = -LjkR,Q,kl2 + 0{h^) 

= -L,,Bj^l2 + 0{h% (73) 

where, because the probe begins at rest, we have taken 
i-''^ = 0(h). Thus (72) generalises this equation to in- 
clude the effect of the quadrupole moment of the parti- 
cle. Because this quantity is time dependent, this allows 
for the possibility of resonance between the probe and 
the wave, the consequence of which we shall explore in 
the following example. 



To argue this more rigorously, suppose that the incident wave 
has amplitude B and frequency Q,, and that the internal motions 
of the probe have frequency ui and amplitude I. We require 
-B 1 in the linear approximation, f2Z <IC 1 to ensure that the 
probe is much smaller than the gravitational wavelength, and 
uil < 1 so that the internal motions are not superluminal. It 
follows from (70) that the first term AX/M ~ BfiVlLO, and from 
(12) that the second term XM/M2 = {X/M){^ d°'Taod^x/M) ~ 
(X/M)BZ2n^2 rpj^g x/M = f AXdt/M can be no larger 

than (AX)inaxAt/M ~ EPVlujAt, where At is the duration of 
the interaction. From this we conclude that the second term 
XM/M^ < B^l^n^ui^At is negligible in comparison with the 
first unless the wave and the probe interact for a very long time 
At ~ {Bl^Cluj'^)~^ . This becomes completely impossible as the 
length-scale of the probe I — > 0. 




• = m/2 



FIG. 1. A toy model detector: two masses, connected by a 
light rod, rotate in the xs = plane; a gravitational plane- 
wave, propagating in the ii-direction, disturbs its centre-of- 
mass in both the longitudinal and transverse directions. 

Example: Rotating Rod 

Let us consider the probe depicted in figure 1, a light 
rod (length 21) with bobs of mass m/2 at each end, spin- 
ning with angular frequency uj about the x^-axis. A valu- 
able feature of equation (72) is that one only needs the 
unperturbed motion of the detector (as captured by lij 
and Lij) to calculate the motion of the centre-of-mass to 
lowest order in /i^^ ; this is not true of a "first principles" 
approach to the problem (see appendix A) which com- 
plicates that calculation considerably. The unperturbed 
locations of the two masses are, in the background, 

x(^) = ;(coswt,sina;t,0) = -x^^\ (74) 

and assuming that the speeds are not relativistic (for the 
sake of simplicity) it is easy to confirm that 

/ - sin(2a;t) cos(2wt) - 1 0\ 
iij — Lij = muP cos(2ajt) -f- 1 sin(2a;i) 

V 0^^. 

Inserting this into (72) and setting the total mass/energy 
M = m under the nonrelativistic assumption, we con- 
clude that centre-of-mass of the spinning rod moves ac- 
cording to 

= '^BlI{t){co8{2ut) + 1), (75) 

in the transverse directions i = 2,3. For a generic grav- 
itational wave, this equation predicts an oscillation in 
the centre-of-mass that averages to zero over many wave- 
lengths. If the wave is of frequency 2w, however, a res- 
onance occurs in which the detector can steadily "walk" 
in the transverse-direction. A gravitational wave of the 
form 

gives rise to an average transverse velocity 
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One of the most surprising aspects of this phenomenon 
is that the walking motion (77) is not associated with 
any transverse momentum: P2 = P3 = 0. The detector 
moves without being pushed, as it were: due to a careful 
conspiracy between the probe's internal motion, and the 
stretching and squeezing of space, the centre of the probe 
is displaced with each period. 

To understand this on an intuitive level, let us imagine 
for a moment that the rod joining the masses does not 
exist, but that at i = the masses have the same posi- 
tions and velocities as before. Because the gravitational 
wave is invariant under translations in the transverse di- 
rections, the transverse momentum (i.e. the transverse 
components of the momentum covectors) of the two par- 
ticles will be conserved, and hence the total transverse 
momentum remains zero. However, the velocity vectors 
of the masses are related to their conserved momentum 
covectors by the physical metric gab, which varies in the 
x^-direction. Thus, because the physical metric differs 
between the positions of the two masses, while their mo- 
menta are equal and opposite, their velocities will not be. 
In this fashion, a gradient in the gravitational field across 
the detector can cause a drift in the centre-of-mass of the 
system. The role of the rod in our detector is simply to 
apply equal and opposite forces to the masses (again, 
having no effect on the total transverse momentum) so 
that once t = tt/2uj and the gradient of the gravitational 
wave across the detector has reversed, the masses are now 
at the same value of x^, and the drift that has occurred 
in the first quarter-wavelength will not be undone. 

In appendix A we substantiate this intuitive picture 
with a detailed rederivation of equation (70) from first 
principles. Not only does this further aid our understand- 
ing of the phenomenon, it should assuage any concerns 
that this unfamiliar effect might simply be an unphysical 
artifact of our formalism. In fact, the subtlety and com- 
plexity of this calculation emphasizes the computational 
advantage of our approach, not only for MoE, but for 
angular momentum and energy-momentum also. 

VI. GRAVITATIONAL PLANE- WAVES 

As a final exploration of our formula (32) for gravita- 
tional intrinsic spin, we shall evaluate s"^^ for a plane- 
wave. The motivation for this endeavour is to point out 
a number of interesting features, and to allow for a com- 
parison with other descriptions of gravitational angular 
momentum. 

A transverse-traccless gravitational plane-wave 

hfj.u ^ hfii/ikaX ), — (1, 1,0,0), 

/i^o = h^,i = h = 0, (78) 

has an extremely simple spin tensor: 

Ks"23 = fc"(/^x/i+-/i+/lx), (79) 



where hx = and h-^. = /i22 = —^33 are the transverse 
components of the wave. As one would expect, s"^^^ de- 
scribes transverse spatial spin flowing in the direction of 
propagation of the wave. Furthermore, the amplitude of 
5^23 quantifies the internal spinning motion of the field, 
as can be seen when wc consider a monochromatic wave 
where the "plus" and "cross" polarisations differ by a 
phase 6: 

= cos {u][t — x^)) , 
hx ^ Ax cos, {uj{t- x^) - 0) . (80) 

In this case, the spin-density is constant over spacetime, 

Ks°23 = ujAxA+sinO, (81) 

and is greatest in magnitude when the wave is circularly 
polarised, that is. when 6 ~ ±7r/2. Note that a wave with 
a purely linear polarisation will carry no spin at all. 

In [1] we saw that the energy-momentum tensor of a 
TT gravitational plane- wave was independent of the time- 
like vector m'^ that defines the wave's TT-frame (8). A 
similar property holds for the spin tensor, but it is com- 
plicated by the fact that spin is constrained to be spatial 
with respect to the TT-frame, that is, m'^s"^j^ = 0. As we 
shall see, the longitudinal and non-spatial spins do trans- 
form as the TT-frame is changed, and in doing so they 
adapt the spin tensor to obey the spatial constraint for 
the new u^; however, the transverse spatial spin current 
s"23 is left invariant. To demonstrate this invariance, we 
perform a gauge transformation on the field (78) that 
maintains its plane-wave form, 

<5V - 5(M mika^x")) = 2A:(^^,), (82) 
and note that the spin tensor changes by 

+ - , (83) 

confirming that 6s"23 = 0- 

Now suppose that the gravitational field (78) has been 
transformed to a new TT-frame, so that in some other 
Lorcntz coordinate system {cc^ } we have /i^'o' = /i^'i' = 
0. Then by the same calculation that led us to (79) 
the transformed spin tensor s'^j^^ will obey s'"^'o' = 
s'" If = exactly as the original tensor did in the 
original coordinate system. The only non-zero compo- 
nent of the transformed tensor (in the primed basis) will 
be s'"2'3', and this quantity will also be gauge-invariant 
by the same argument we used for 5^23 • These two 
gauge-invariant currents arc related by the constant fac- 
tor 2A[^2A'^l3, where A^i, is the Lorentz transformation 
between the two coordinate bases: 

a 'a \ n' a ly' la. 

S 23 = S 23 = ^ 2A 3S 

= (2A[2'2A3'l3) s'%,3,. (84) 
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This constant of proportionality ensures that s°j^ and 
■*'"bc describe exactly the same spatial transverse spin 
current in either basis: 5^23 = ^'"23 ^^'^ = s'^2'3'- 

Thus, the only effect of a change in TT-frame is to re- 
express the same physical information (the transverse 
spin current of the wave) in terms of spin that is spa- 
tial with respect to a new rest-frame. In the absence of 
some material body (a detector or a source, for exam- 
ple) the massless gravitational plane-wave cannot define 
a preferred rest-frame, and so the spatial nature of its 
intrinsic spin will always have this ambiguity. 

As a consequence of this, while a plane- wave region can 
"sew together" two different TT-frames to form a seam- 
less picture of the propagation of gravitational energy- 
momentum (as described in section III D of [1]) the same 
cannot be done for angular momentum: there will always 
be a discontinuity where the spatial spin of one frame is 
converted into the spatial spin of the other. Even so, one 
can construct a gravitational spin pseudovector 

= e"^''''sAM./2, (85) 

which is truly independent of TT-frame, and will there- 
fore give a continuous description of gravitational spin 
within the sewing region. The invariance of s" follows 
directly from the totally antisymmetric part of (83): 
Ss[ap,v] = 0. The physical interpretation of this pseu- 
dovector is not immediately clear, but suffice it to say 
that for a plane-wave, s" captures only the spin that is 
linearly independent of the wave- vector fc'^.^^ For the 
plane-wave (78) we have been studying, the spin pseu- 
dovector is 

Ks" = k"{h^h+~h+h^), (86) 

capturing all the physically pertinent information of (79) 
in a completely frame-independent fashion. 

Finally, we should highlight the major difference that 
exists between the gravitational spin currents in (79) and 
the corresponding quantities given by the traditional ap- 
proaches, including the Landau-Lifshitz tensor [4] and 
the integrand of the ADM energy- momentum [5]. In 
these descriptions, the local energy-momentum and spin 
of the gravitational field are packaged together in a sin- 
gle object, a Belinfante energy- momentum tensor t^,^ ^ 
dhdh + hd^h + 0{h^).^^ The local angular-momentum 
currents are then x^^t^]" alone, with no extra "intrin- 
sic" component. According to this viewpoint, there is 



We also note that s" bears a resemblance to the Pauli-Lubanski 
pseudovector S" = e"^'"'P;^L^y/2, which characterizes the total 
spin of a particle or matter field, and reduces in the particle's 
rest-frame to (mass times) the familiar axial angular-momentum 
vector of non-relativistic mechanics [18] . 

A Belinfante energy-momentum tensor can be constructed from 
any energy-momentum tensor and spin tensor, including our 
own: tfj,^[T,s] = + da(Sf_,„°' + s^ij" - s°'^^)/2. We perform 
this calculation in [2] and compare the result with the Landau- 
Lifshitz and ADM Belinfante tensors discussed here. 



no transverse angular momentum within a harmonic- 
gauge plane- wave: X[2t3]" = 0.^^ This differs dramat- 
ically from our description (79) and stands opposed to 
the intuitive notion of intrinsic spin as quantifying the 
internal spinning motion of the field. Without separating 
gravitational energy-momentum and spin into two sepa- 
rate tensors, t^^, and s"^j^, the intrinsic spin carried by 
a (harmonic-gauge) plane-wave can never be manifestly 
present within the wave. 

To be clear: the Belinfante-style descriptions still cor- 
rectly quantify the total angular momentum of the wave, 
but they assign this angular momentum to the wave's 
boundary, not its interior. Considering that the angu- 
lar momentum currents along this boundary are given by 
^[2^3]" as always, and are thus explicitly dependent on 
x'^ , even these currents cannot be thought of as a local 
and intrinsic property of the field. This perverse picture, 
in which all the spin of a gravitational wave resides on 
the edge of the wave, and this supposedly intrinsic quan- 
tity depends on the coordinate distance from the ori- 
gin, only emphasizes what was already well-known: the 
Landau-Lifshitz tensor and the integrand of the ADM 
energy-momentum should not be taken seriously as local 
descriptions of gravitational energy-momentum or spin. 
While they certainly define meaningful global quantities 
[19], the gauge-freedom of these Belinfante tensors cannot 
be fixed in a natural manner, and they commonly display 
negative energy-density and spacelike energy-flux. 

VII. CONCLUSION 

Together, the energy-momentum tensor T^,y and the 
spin tensor s"^^ completely characterize the energy, mo- 
mentum, and angular momentum carried locally by the 
linearized gravitational field: 

= jd^h^pd.h''^ , (87) 

^s'^^,^2hpi,dn^]^l (88) 

The gauge freedom of this description is highly con- 
strained by the harmonic gauge condition, 

V = o> (89) 



This follows from simple index combinatorics. Within the plane- 
wave t/j,^ ~ kkhh + kkhh + 0{h^), and because k^k^^ = and 
k^h/^i, = 0, both the free indices must occur on the wave- vectors, 
i.e. tfii, oc kf^kjj. This continues to be true at higher order, where 
the terms in t^^ are of the form kkhhh"-' + kkhh"-^ . Con- 
sequently, the transverse angular momentum vanishes exactly: 

To avoid a discussion of the boundary at infinity, suppose the 
plane-wave is in fact restricted to a spatially compact region; 
in this case, one will find that 3;[2i3]°' ^ at the boundary of 
the region, and the spatial integral of aipta]'' will amount to the 
same total angular momentum described by s''23- fact, it is 
generally true that (under suitable boundary conditions) t^i, [r, s] 
gives the same global measure of energy-momentum and angular 
momentum as T^,y and see [2] for details. 
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which arose as a consequence of the derivation of t^^ ; the 
last remnant of this freedom is then ehminated by insist- 
ing that the incident gravitational field be transverse- 
traceless, a program motivated in part by appealing to 
the gauge-invariant exchange of energy-momentum be- 
tween gravity and an infinitesimal probe, and also dis- 
tinguished by the numerous desirable properties that the 
tensors display in transverse-traceless gauge: positive 
energy-density, causal energy-flux, and spatial spin. 

We developed this framework around a simple princi- 
ple: wherever the energy, momentum, or angular momen- 
tum of matter is changed under the influence of gravity, 
there must be an equal and opposite change in the energy, 
momentum, or angular momentum of the gravitational 
field. This idea, and the requirement that r^^ be sym- 
metric and free of second derivatives, was enough to de- 
termine the energy-momentum tensor (87) and the field 
condition (89). To determine the spin tensor uniquely, we 
demanded that it satisfy two physically-motivated con- 
ditions: first, the spin of a gravitational plane-wave must 
flow in the direction of propagation of the wave (24); sec- 
ond, a transverse-traceless field must possess a traceless 
spin tensor (31). The latter condition ensures that local 
field redefinitions suffice to cast t^^, and s"^^ as sources of 
gravity in a quadratic approximation to general relativity 
(29) and simultaneously rids the gravitational field of in- 
finite pressure gradients (see appendix B). The resulting 
spin tensor (89) displays a number of notable properties 
that were not required of it: the Newtonian potential has 
vanishing spin-tensor (28) and any transverse-traceless 
field carries purely spatial spin (27). 

The microaverage, which defines the gauge-invariant 
exchange of energy-momentum between gravity and an 
infinitesimal probe, also renders the exchange of spa- 
tial angular momentum gauge-invariant (47) without the 
need for gauge-fixing. In the same system, the ex- 
change of non-spatial angular momentum can displace 
the center-of-mass of the detector, beyond that which 
would be expected due to the exchange of linear mo- 
mentum alone (72). Indeed, if the internal motions of 
the probe resonate with the incident wave, the detector 
may "walk" in a transverse direction, and acquire a net 
displacement over many wavelengths. We have explored 
this phenomenon for the specific example of a rotating 
rod (75) and rederived our predictions from first princi- 
ples (see appendix A). 

Unlike r^^, the spin tensor of a gravitational plane- 
wave is not completely independent of the TT-frame in 
which the wave is prepared. While the current-density 
of transverse spatial spin (in any frame) is invariant, the 
full tensor adapts so as to remain spatial with respect to 
whichever TT-frame is used. Thus, if a plane-wave region 
is used to sew together two TT-frames and produce a 
seamless picture of energy-momentum propagation, there 
will inevitably be a discontinuity in s"^^, where the spin 
is projected from one spatial hypersurface to another; 
however, a spin psuedovector can be defined (85) that is 
conserved across this interface. 



The spin carried by a plane-wave (79) is also an excel- 
lent example with which to compare our framework to 
the familiar "Belinfante" energy-momentum tensors of 
Landau and Lifshitz, and Arnowitt, Deser and Misner. 
Whereas s"^^ describes spin that is manifestly present 
within the wave, the density of which depends on the 
rotational motion of at each point, the Belinfante 
tensors assign all angular-momentum to the boundary of 
the wave, and its density there is not simply a function of 
h^i, (as a truly local intrinsic property of the field would 
be) but is also dependent on the distance of the point 
from the origin. 

Returning to our previous paper [1], it becomes clear 
that many of the remarkable properties of our gravita- 
tional energy-momentum tensor (including its positive 
energy-density and causal energy-flux) owe their exis- 
tence to the careful separation of gravitational energy- 
momentum and gravitational spin. Now that we have 
made this separation explicit, and derived a formula for 
s"^y , we have all the ingredients necessary to understand 
the broader theoretical picture in which our description 
resides. This is the task of our next paper [2], the results 
of which, in many respects, are the main reward for our 
work here. 
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Appendix A: Moment of Energy Exchange from 
First Principles 

In order to rederive equation (70) from first principles, 
we shall consider a detector, centred at the origin, com- 
posed of a set of N test-particles connected by some form 
of "light" mechanical apparatus. "^^ The n"^ particle has 
mass m„, proper time r„, and follows a worldline 2/^^(t„) 
in the physical spacetime; its 4- velocity = dj/^/dr„ 
has unit norm: u'^u^gab = "1- In this approach, the de- 
tector is not truly infinitesimal, but we stipulate that the 
length-scale of the detector Z ~ j/^ be sufficiently small 
that we may ignore terms 0{l^) in our calculation, leav- 
ing us with a quadrupole approximation of the probe. As 
usual, a weak gravitational plane-wave is incident upon 
the detector, represented in transverse-traceless gauge 
in the background: /i^^ = hj^^™ as given in (62). As 



The adjective "light" is used to indicate that the total energy- 
momentum of the apparatus is negligible compared to that of 
the particles. 
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d2h^i, = d^h^i, = 0; the physical spacetime is isomet- 
ric under translations in the and directions of the 
{y^} coordinate system; for this section, we will reserve 
the index fc = 2, 3 for these transverse directions. 

In the physical spacetime, the energy-momentum ten- 
sor of the particles is 



(-^particles) b 




(Al) 



where pna = m-ngabu'^ is the 4-momentum of the n'^ par- 
ticle. For the purposes of defining the moment-of-energy 
of the detector, we assume that the energy-momentum 
of the light apparatus is negligible: 




particles 



(A2) 



In terms of background quantities, this is 

N 



Xi= - ^ m„(770Q -I- Koa)x'"Xr, 



n=l 
N 

EiO 

n=l 



(A3) 



where x^(t„) are the coordinates of the particles in 
the background spacetime, primes indicate differentia- 
tion with respect to t„, and = /ip^(a;^(r„)) is the 
gravitational field evaluated at the n*'^ particle. The rate 
of change of the moment-of-energy is therefore 



N 



Xi = ^ m„ 



71 = 1 



The normalisation of the 4- velocity. 



-1 



'ix„ ) + {hnij + ^ij)XnXn 1 



(A4) 



(A5) 



ensures that x'° ~ 0{1) and x'^° ~ 0{P), so the first 
term (A4) is 0{P) and can therefore be neglected. Con- 
sequently, 



X,. 



N „ 
™/0 



Oil' 



(A6) 



The accelerations a;"j in (A6) arc caused by both the 
gravitational field and the mechanical forces exerted on 
the particles by the apparatus. Our aim is to infer Xi 
without assuming any detailed model of the latter. This 
might seem an impossible task, as it appears that we will 
need to know the motions of the particles (or the forces 
from the apparatus) to first order in h^,^ if we wish to 



calculate the first order contribution to Xk- Fortunately, 
because the apparatus is light, and the transverse mo- 
mentum is conserved, only the unperturbed motions of 
the particles will be required. To see this, we start by 
calculating the linear momentum of the probe: 




.g(T'particlos)°id^2/ 



(A7) 



n=l 



where once again we assume that the momentum of the 
apparatus can be neglected. Because the physical space- 
time is isometric under translations in the and y' 
directions, the transverse momentum will be con- 
served:'^^ 

= Pk=dt m„(4i + hnkt)x'^^ 

N ^ 



(A8) 



which is equivalent to the statement that the mechanical 
forces on the particles (due to the apparatus) balance 
one another. ^'^ Substituting this constraint into equation 
(A6) we find that 



N 



Xk = -J2^drJKk^x':^ + Oif) 
,1=1 

= dt (^-J2^mnhnk^x'^^ +0(1^), 



(A9) 



which is easy to integrate 

N 



34 



Xk 



(AlO) 



This follows from the standard argument: dhSap = guarantees 
that {dkT is a Killing vector, V^'^{dkf^ = 0, and thus = 
\f-9^ a{T\(dk)^) = da{\^T\), the spatial integral of which 
is Pfe = 0. 

Although the total momentum of the apparatus is assumed to 
be negligible, we have not made any assumptions about the local 
flux of momentum between the apparatus and the particles, and 
so the individual mechanical forces on each particle cannot be 
neglected. The constraint (A8) arises because the apparatus has 
much less mass than the particles, and so any momentum it were 
to gain would send it off with a very large velocity that would 
be impossible to maintain while in contact with the particles; in 
order to stay connected to the particles, the momentum of the 
apparatus must remain very small, and the forces acting on the 
apparatus must (approximately) sum to zero. 
■^■^ The constant of integration is set to zero by the initial conditions: 
the probe is at rest {Xi = 0) before the wave arrives (/x^i^ = 0). 
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This is the equation we sought: every instance of x'^ is 
multiphed by h^^,, so only the unperturbed motions are 
needed to determine to linear order in the gravita- 
tional field. 

The last step is to relate the hnki to the gravitational 
field at the origin: 

hnt^u = hf,,^{t, 0) + X'„^^hf,^.{t, 0) + 0{l^) 

= Bll{t)-xli3ll{t)+0{n (All) 
as a result, equation (AlO) becomes 

= -BY, m„<:^ + i3Y + 0{l^) 

= -BY {Xk) +BY (j2^m,^<^':^ +0(1') 
= BY m„:.i,<^ + Oih') + 0{P). (A12) 
This simplifies even further when we notice that 
hj - =dt[- I V^iTpaiticicsfoyiyjd^y 

2 / \/^(rparticlos)°[jyj]d^y 




N 



dt i^nin 

\ti=1 

+ 0{h) 



2E 



n=l 



N 



N 



n '^n 



2 ^ ^ TTljixJ^ X'^ 2 ^ ^ TTLjiX^^ 
n—1 n—1 

+ 0{h) + 0{l^) 

N 

2^m„x^„< + 0(/i) + 0(/*), (A13) 



which gives us our final result: 

Xk = BY - La) /2 + Oih') + 0(/3), (A14) 

exactly as predicted by equation (70). 

It should be clear that our formalism provides a much 
more direct route to this result: one needs only to pro- 
duce (64) and integrate, a straight-forward operation 
that lacks the "insightful" steps of the first principles 
calculation, such as invoking conservation of transverse 
momentum (A8) to remove degrees of freedom from (A6). 
However, the moral of this appendix is not simply that 
our method is more computationally efficient; equally im- 
portant is the intuitive power that our framework con- 
fers. Working from first principles, it is hard to imag- 
ine that one would have thought to derive (A14) in the 
first place, as there is no obvious reason to expect that a 
gravitational wave would produce a transverse motion in 




FIG. 2. The torques on an infinitesimal cube of vacuum due 
to the flux of gravitational intrinsic spin. 



the detector's centre-of-mass. In comparison, our unified 
picture of local gravitational energetics brought this phe- 
nomenon to mind as readily as the exchange of energy, 
momentum, or angular momentum. 



Appendix B: Physical Interpretation of the 
Traceless Condition 



Here we offer an alternative derivation for the traceless 
condition (31) which helps to illuminate its physical inter- 
pretation. Let us begin by examining the spin tensor of 
the transverse-traceless gravitational field in detail. The 
purpose of this analysis is to isolate an algebraic prop- 
erty of s°^j, that signifies unphysical behaviour, and then 
design a condition so that this possibility cannot arise. '^^ 

Because the spin of a transverse-traceless field is spatial 
(Property la) we can write 



a a 



(Bl) 



where s"^ is the axial spin tensor, the current-density 
of intrinsic spin about the a;*-axis. Each component 
represents the fiux of s^-axis spin in the x^-direction; in 
other words, SyS is the torque (along the a;-' -axis) that 
acts on a small surface (x' = const.) of area S. 

Let us consider the I — > limit of an / x Z x ^ cube 
of vacuum {T^i, = 0, hfj,i, ^ 0) as depicted in figure 2. 
The torque along the cc^-axis, acting on the = face, 
is G\ = S12P ■ There will also be contributions from the 
X[iTj]k part of the angular-momentum current density, 
but these terms will be of order and so can be safely 



Note that we restrict our attention to the spin of the transverse- 
traceless gravitational field. The spin tensor can only be ex- 
pected to have a sensible physical interpretation under the same 
conditions that r^^ describes positive energy-density and causal 
energy-fiux, i.e. for all TT-fields, arbitrary (harmonic-gauge) 
plane-waves, and static fields. It will be trivial to extend Condi- 
tion 2 to include the last two cases, and since their inclusion does 
not constrain A\ , it is simpler to ignore them in what follows. 
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neglected. It is convenient to think of Gi as being gener- 
ated by two equal and opposite forces Fi = 2si2l acting 
on the points (0, Z/2,Z/4) and (0,^/2, 3Z/4) as shown in 
the diagram. On the opposite face (x^ ~ I) there will be 
a torque along the x^-axis G'l = — (si2 + ldiSi2)P, the 
minus sign arising as a result of the opposite direction 
of the outward normal, and the second term being neg- 
ligible as long smooth in the cube. Again, this 
torque can be thought of as being generated by equal 
and opposite forces Fi acting at the points (/,//2,//4) 
and {1,1/2,31/4:). Following the same approach, we ren- 
der the x^-axis torques on the = and x^ = I faces 
as forces F2 = 2s2il acting on the appropriate points on 
the cube. 

We now split the cube along the plane x^ = 1/2, and 
consider the two "half-cubes" separately. The isotropic 
pressure acting on each half-cube can be evaluated using 
the formula 



P 



6 



1 ^ 



fr ' 



(B2) 



where the index r enumerates the six faces of the half- 
cube (each with area Ar and outward unit normal n^) 
and fr is the force acting on the r"^ face."^^ For the 



upper half-cube {x^ > 1/2) both Fi forces are directed 
inwards, while the two F2 forces are outwardly directed; 
thus (B2) gives 



-I2F2-2F1 4(S12-S2l) 



6 



P/2 



31 



(B3) 



where we have once again ignored the negligible forces, 
such as TasZ^ on the x^ = I and x^ — 1/2 faces. The 
calculation for the lower half-cube {x^ < 1/2) is identical 
except that the forces Fi point outward and F2 point 
inward; as a result, Piowcr = — ^uppcr- Therefore, within 
the cube we find a pressure gradient 



dP 

dx^ 



Pu 



1/2 



16(si2 - S21) 
3/2 



(B4) 



which grows without bound as the limit Z — > is taken! 
The only way to avoid these infinite pressure gradients is 



to insist that s 



[jj] = 0, or equivalently 



0. 



(B5) 



Hence the traceless condition (31) can be derived by re- 
quiring that the (gauge-fixed) gravitational field be free 
of infinite pressure gradients. 
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